function y=sigmoid(x)

y=1./(1+exp(-x));

return